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ABSTRACT 

We present a technique for constructing equilibrium triaxial TV-body haloes with nearly ar- 
bitrary density profiles, axial ratios and spin parameters. The method is based on the way 
in which structures form in hierarchical cosmological simulations, where prolate and oblate 
haloes form via mergers with low and high angular momentum, respectively. We show that 
major mergers between equilibrium spherical cuspy haloes produce similarly cuspy triaxial 
remnants and higher angular-momentum mergers produce systems with lower concentrations. 
Triaxial haloes orbiting within deeper potentials become more spherical and their velocity dis- 
persion tensors more isotropic. The rate of mass loss depends sensitively on the halo shape: a 
prolate halo can lose mass at a rate several times higher than an isotropic spherical halo with 
the same density profile. Subhaloes within cosmological simulations are significantly rounder 
than field haloes with axial ratios that are ~ 30% larger. 

Key words: methods: iV-body simulations - galaxies: haloes - galaxies: interactions - cos- 
mology: theory - dark matter. 



1 INTRODUCTION 

A generic prediction of the currently favoured cold dark matter 
(CDM) cosmological model of hierarchical structure forma- 
tion is that dark matter (DM) haloes of gala xies and c lusters 
are flattened triaxial systems (e.g., iBarnes & Efstathioil 11987 * 



i Frenk et alJll988tlDubinski & Carlberg|'l99 lHWarren et al]ll992t 
Icole & Lacevl Il996t 'Thomas et alJ 1 1998k Ijing & Sutcl hool) . 

Observationally, inferring the intrinsic shape of DM haloes 
of galaxies and clusters is a difficult task, but nonetheless 
several potentially powerful probes exist that may allow us 
to to distinguish between spherical and flattened DM haloes. 
These include the dynamical modelling of collisio nless tracers 
such a s tidal stre a ms or bi ting the Milky Wa y Jjohnstoii^yilJ 
1999[ llbata et alJ 1200 ll iMaver et alJ bool iMaiewski et all 
2004), the distri butio n and k ine matics of gas in spiral ga laxies 

ikeji^^r^maineJ__[^94^__[£ranxjanGoA 
1994 ; Ischoelnri^erZFraiix&d^ ^97^ [Merifiek 



2002), and polar ring galaxi es IschwejzeZ^hitoo^^^^ir 
198a 'Sackett&SparkeNl99d: ISparkd 12001 llodice et alJ 120031) 



which provide shape constraints perpendic ular to the disc 
plane, gravitational lensing applications I Kochaneja P995 ; 
B^rtejmami ] _StejrmTejz & Weissl Il995t liCoopmans et alJ 1998 ; 



Qguri^ee^^Sutoll2003ir and the flattenin g of the extended X-ra' 



isophotes in elliptical galaxies 'Buote & Canizaresll 19941 1 19961 
1998: lBuoteetall2002l) . 

Numerical experiments of isolated equilibrium models are 
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very useful for studying the dynamical evolution of gravitat- 
ing systems in a controlled way. Several techniques exist for 
constructing iV-body realizations of spherical haloes and multi- 
component galaxies iHem quistj |l993|: iHollev-Boc kelmann et alJ 
l2003t iKazantzidis. Magorrian & Moorell2004al hereafter KMM), 
but it is far more complicated to build triaxial equilibria. Even 
though the Jeans theorem guarantees that the distribution function 
depends only on the isolating integrals of motion, explicit expres- 
sions for the latter other than the energy per unit mass, E, are 
rarely known in the case of triaxial potentials. As a result, the tri- 
axial models that have been constructed so far have been either 
limited to a few special analytical cases (e.g., Stackel potentials or 
rotating f{Ej) models with E.j being the Jacobi constant) or en- 
tirely based on numerical techniqu ej_jSchwamcMdJ^^79j, [T993t 
iMerritt & Fridmarll996tlPoon & Merrittl200 j:lTerzicl2003t) . 

Recently, iHollev -Bockelma nn et al] 1200 ll) used a technique 
of adiabatically applying a drag to the velocities of the parti- 
cles along each principal axis in order to cr eate cuspy triax- 
ial systems starting with a spherica l Hernquis t |Hem auistlll99ol) 
model. T he original prescription o fjHernc[uisjl E^l) was gener ~ 
alized by Boily. Kroupa & Penarrubia-Garrido (2001) for accomo- 
da ting composite, axisymr netric models of galaxies and extended 
bv lTinker & Rvdei] {2002 ) for investigating the effect of rotating, 
triaxial halos on disk galaxies. These techniques are useful and have 
been used to understand the response of systems to live triaxial po- 
tentials. However, they are somewhat restricted to modest values of 
the flattening and it is difficult to incorporate a significant amount 
of rotational flattening. 

Triaxiality can arise in a number of ways. The most common 
way of constructing numerical models of triaxial galaxies is via 
the merger of other objects. Examples include binary mergers of 
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spherical haloe s (e.g.JWhitdl978llFurton & BarnesfeOOlh and disk 
alax i es ("e.g.. lGerhardll98lt |Barneslll992t iBarnes & Hemauistl 



1 + c 



3c r vir 



(3) 



Il996t lNaab & Burkerll l2003t lkazantzidiset alJl2004bh. as well 
as multiple me rgers of systems (e.g., |Weil & Her nauistl Il996t 
lDubinskilll998h . The structure of the final remnant of two com- 
ponent models depends sensitively on both the orbital geometry 
<Naab & Burkertll2003l) as well as on the inclination and internal 
properties of the disks iKazantzidis et alj20 04b). 

A general technique for constructing cuspy axisymmetric and 
triaxial TV-body systems would be advantageous for many purposes 
including studying the dynamical friction in flattened and rotating 
systems, the mass loss and tidal stripping from triaxial substruc- 
ture haloes, the effects of baryonic accretion and disc formation on 
triaxial halo shape and anisotropy, the effects of triaxial shape on 
the properties and stability of discs, the interaction between central 
black holes and cusps and the inflow of gas in triaxial systems. 

In this paper we explore the generation of triaxial structures 
that is based on the way that haloes in hierarchical models of struc- 
ture formation obtain their shapes. Inspection of a cosmological 
simulation reveals that triaxiality arises via mergers that take place 
with various amounts of angular momentum which is generated 
from the large-scale tidal field. Mergers between two haloes that oc- 
cur with little angular momentum (radial mergers) produce prolate 
systems, whilst high angular momentum mergers produce oblate 
systems. Most CDM haloes form via a sequence of mergers with 
varying amounts of angular momentum such that haloes with ar- 
bitrary triaxiality may be formed and the triaxi ality will vary with 
radius I M oore et all200lt IVitvitska et all2002l) . 

The outline of this paper is as follows. In Section 2, we de- 
scribe our technique for constructing equilibrium cuspy TV-body 
haloes with various degrees of flattening. There we discuss in detail 
the numerical experiments we performed and present our results for 
the internal structure of the resulting models. In Section 3, we in- 
vestigate the tidal evolution of triaxial substructure haloes within 
a static host potential and study the shape of subhaloes within a 
cosmological CDM simulation. Finally, we summarise our main 
results in Section 4. 



2 CONSTRUCTING CUSPY TRIAXIAL DARK MATTER 
HALOES 

We begin with a cuspy sp herical DM halo which follows the 
iNavarro. Frenk & White! fl996. hereafter NFW) density profile 



p(r) = 



(r/r s )(l + r/r a ) 2 



(r < r vir ). 



(1) 



Here the characteristic inner density p s and scale radius r s , are sen- 
sitive to the epoch of halo formation and tightly correlated with 
the halo virial parameters, via the concentration, c, and the virial 
overdensity A v i r . Since the NFW density profile corresponds to a 
cumulative mass distribution that diverges as r — > oo, we impli- 
ment an exponential cut-off for r > r rv i r which sets in at the virial 
radius and turns off the profile on a scale rd CC ay which is a free 
parameter and controls the sharpness of the transition 
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where c = r v i r /r s is the concentration parameter. Finally, in order 
to ensure a smooth transition between Q and at r v ir, we require 
the logarithmic slope there to be continuous. This implies 



Once the density profile p(r) of our initial model has been 
specified, we construct an TV-body realization using the technique 
discussed in detail in KMM. To summarize, the particle positions 
and velocities are initialized by sampling the exact phas e-space dis- 
tribut ion function using Eddington's inversion formula (Eddinaton 
1916) for the given density profile p(r) and assuming spherical 
symmetry. Note that in this paper we will consider only initial mod- 
els with isotropic velocity dispersion tensors and therefore their 
distribution function will only depend on the binding energy E, 
f — f(E). For the following experiments we begin with an ini- 
tial galaxy-sized halo with virial mass M V ir = 7 x 10 h~ Mq 
and concentration of c = 10 in accordance with theoretical pre- 
dictions for objects at this mass scale and for the adopted cosmol- 
ogy (e.g.. iBullock et af]|200ll) . Note that for all simulations pre- 
sented in this paper including the cosmological ones described in 
the next section, we consider as our framework the concordance 
flat ACDM model with present-day matter and vacuum densities 
fi m = 0.3 and Qa = 0.7, respectively, dimensionless Hubble con- 
stant h = 0.7, present-day fluctuation amplitude as = 0.9 and 
index of the primordial power spectrum n — 1.0. For this choice of 
concentration, the virial and scale radius of the initial halo model 
are r v i r = 181 /i -1 kpc and r B = 18.1 h -1 kpc, respectively. The 
initial model is constructed with TV = 5 x 10 5 particles and we use 
a gravitational softening of e = 0.3 hT 1 kpc. 

All the simulations presented in this paper were carried out 
using PKDGRAV, a multi-stepping, parallel TV-body code (Stadel 
2001) which uses a spline kernel softening and multi-stepping 
based on the local acceleration of particles. The time integration 
was performed with high enough accuracy such that the total en- 
ergy was conserved to better than 0.1% in all cases. We have also 
explicitly checked that our results are not compromised by choices 
of force softening, time-stepping or opening angle criteria in the 
treecode. 

We perform a first experiment in which we merge two iden- 
tical spherical haloes with zero angular momentum corresponding 
to a head-on collision with an impact parameter 6 = 0. The ini- 
tial relative separation of the two haloes is r scp = 210/i -1 kpc 
and are they infalling at a relative speed of v TC \ = 300 kms -1 . 
Relative velocities significantly higher than this value lead to un- 
bound orbits, so we chose this value so as to maximize the amount 
of flattening from a single merger. Note that the initial separation 
corresponds to haloes which are overlapping. This does not influ- 
ence the final result as nearly identical remnants were produced 
with non-overlapping initial models. The two haloes merge within 
a time-scale of order the crossing time t C ro SS = \J r^JGMvii 
which is approximately equal to 2 Gyr for our initial models. In all 
simulations the merger remnants were allowed to settle into equi- 
librium for several crossing times after the merger was complete 
which was established by monitoring the fluctuations of the den- 
sity profile of the final system. In addition, the virial ratio 2T/W 
varied between 0. 1 and 0.3% from 1.0 for the same time-scales con- 
firming that our triaxial remnants are indeed in equilibrium. In all 
cases the centre of the remnant is identified using the most bound 
particle, which agrees very well with the centre of mass recursively 
calculated using smaller spherical regions. 

The radial merger produces a prolate axisymmetric system 
with axial ratios a:b:c = 2:1:1 where a, b and c are the long, in- 
termediate and short axis, respectively. For each remnant, we cal- 
culate principle axis ratios s = b/a and q — c/a (a > b > c), 
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Figure 1. Density profiles of the initial model (dotted curve), the a:b:c = 2:1:1 prolate model (solid curve) and the a:b:c = 2:2:1 oblate model (dashed curve). 
The density was calculated using spherically averaged (left panel) and ellipsoidally averaged (right panel) bins such that the density is projected onto the 
symmetry axis. 
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Figure 2. Circular velocity curves scaled to the same peak circular velocity for the initial model (dotted curve), prolate a.b.c = 2:1:1 halo (solid curve) and 
oblate a:b:c = 2:2:1 halo (dashed curve). The mass distribution was calculated using spherically averaged (left panel) and ellipsoidally averaged (right panel) 
bins. In each case the circular velocities are scaled to the same peak velocity assuming V c oc r. 



from the eigenvalues of a m odified dimensionless inertia tensor 
(e.g.,|bubinski & CarlberJl99lh : ly = ^2 X i x j/<1 2 > where q 2 = 
x 2 + (y/s) 2 + {z/q) 2 is the ellipsoidal coordinate. We use an iter- 
ative algorithm starting with a spherical configuration (a — b — c), 
and use the results of the previous iteration to define the principle 
axes of the next iteration. The procedure is iterated until the val- 
ues of both axial ratios b/a and c/a have a percentage change of 
less than 10 -3 . The dimensionless spin parameter of this remnant, 
which is defined by 



GM2 

where J, E and M are the total angular momentum, binding energy 
and mass respectively, corresponds to A = as expected in a radial 
merger. 

In a second experiment we start with the two spherical haloes 
separated again by r scp = 210/i _1 kpc, but we give one halo a 
transverse velocity equal to the circular velocity of the combined 
model at that distance, u rc i = 240 kms . The circular orbit 
merger produces an oblate halo with axial ratios a:b:c = 2:2:1. As 
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Figure 3. Logarithmically spaced surface density contours of the oblate model viewed in the plane of the long axis (left) and short axis (right). 



we demonstrate later (see Figure [4}, haloes with any amount of 
triaxiality between these values could be created by using unequal- 
mass mergers or mergers with less or more angular momentum. The 
spin parameter of the final system is A = 0.1. As expected, with 
these experiments we have been able to reproduce most of the range 
in A seen in cosmological iV-body simulations (e.g.. Warr en et alJ 
Il992t Lemson & Kauffmann 1999). Although the formation his- 
tories of halos with a hierarchical universe is significantly more 
complex than these binary mergers, we might expect that oblate 
halos have a higher spin parameter and more angular momentum 
than prolate halos. Figure Q shows the density profiles of the ini- 
tial spherical model and that of the prolate and oblate haloes. The 
panel on the left shows results using spherically-averaged bins. The 
density profiles on the right panel were calculated as a function of 
the ellipsoidal coordinate q, where q 2 = x 2 + (y/b) 2 + (z/c) 2 . In 
agreement with other studies (e.g.,|Whitdll978tlVillumsedll983t 



lBarnesll999ljFulton & BarnesfcOOltlBovlan-Kolchin & M^Ooll) . 
remnant central density slopes remain unchanged. In Figure [2] we 
show the circular velocity curves of the initial model and the fi- 
nal prolate and oblate haloes using again spherically averaged (left 
panel) and ellipsoidally averaged (right panel) bins. These have all 
been scaled by the same factor in radius and velocity such that all 
the curves have the same peak circular velocity. The prolate halo 
has a nearly identical concentration to the initial spherical model, 
but the oblate halo is nearly a factor of two less concentrated. 

In order to create haloes with higher amounts of flattening it 
is necessary to merge together the products of the first set of sim- 
ulations. Merging two prolate haloes along the long axis again at a 
velocity such that the system is just bound produces a a:b:c = 3:1:1 
remnant, while merging two oblate haloes in the plane of the long 
axis, but counter-rotating, results in a further flattened system with 
axial ratios a:b:c = 3:3:1 and no net angular momentum, but with 
particles streaming in opposite directions. Figure|3|shows two con- 
toured images of the oblate halo viewed projected along the short 
and long axis, respectively. 

Finally, we performed a minor merger between two spher- 
ical halos with mass ratio equal to 4 : 1. The most massive 
halo was identical to the ones used in the previous experiments 
while the conc entrat ion of the smaller one was scaled according to 
Bullo ck et alJ <200ll) for the adopted ACDM cosmology. The or- 
bital angular moment um of the merger wa s typical of cosmolog- 
ical encounters I Khoc hfar & Burkerlll2003l) . The number of par- 
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Figure 4. Axial ratios and triaxiality parameter of the minor merger remant. 
The solid, dotted and dashed lines, respectively, show results for b/a, c/a 
and T. The profile of the triaxiality parameter shows that this remnant is 
strongly triaxial with T ~ 0.5. The error bars represent the rms variation in 
the measurements from 20 time intervals spanning 10 half-mass dynamical 
times of the remnant at the end of the simulation. The errors in the axial 
ratios and triaxiality parameter are small demonstrating that the remnant is 
in equilibrium. We define the virial radius of the remnant, r v i r , to be the 
distance from the centre where the mean enclosed density is ~ 337 times 
the mean universal density. 



tides and gravitational softening used for the less massive halo 
was such that both halos were resolved with the same mass and 
force resolution. This is useful in the intepreta tion of th e results 
regarding the structure of the remnant. Kazantz ldis et alJ f2004b) 
examined the evolution of disk galaxy merger remants and found 
that their shapes evolve in the outer regions for several crossing 
times. In Figure |4| we demonstrate the stability of this remnant 
by running the simulation for a further 10 half-mass dynamical 
times and calculated profiles of axial ratios and triaxiality parame- 
ter [T = (a 2 — b 2 )/ {a 2 — c 2 )] at 20 intermediate steps. The error 
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Figure 5. The rate of tidal mass loss over a period of 5 Gyr (~ 2 orbital pe- 
riods) for the initial spherical model (dotted line), the a:b.c = 3:1:1 prolate 
model (solid curve) and the a.b.c = 3:3:1 oblate model (dashed curve). The 
prolate halo experiences much more efficient tidal stripping than both the 
spherical and oblate halo on the same external tidal field and orbit which is 
probably caused by the orbital distribution of particles supporting its shape. 
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Figure 7. The anisotropy parameter /3 as a function of radius for the pro- 
late halo model (solid line), oblate halo (dashed line) and initial condi- 
tions (isotropic spherical halo, dotted line). The open symbols show the 
anisotropy parameter for the prolate halo after orbiting for 5 Gyr within a 
deeper potential. 



3 TIDAL EVOLUTION OF TRIAXIAL SUBHALOES 
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Figure 6. The axial ratios as a function of radius for the prolate (solid lines) 
and oblate (dashed lines) 3:1 haloes. The open circles show the axial ratios 
of the same prolate model after orbiting for 5 Gyr within a deeper potential. 



bars in Figure[4]show the rms variation in the profiles at each radius 
and demonstrate that the remnant is in equilibrium. 



Numerical investigations of the tidal stripping and mass loss of 
substructure haloes in CDM models have previously been re- 
stricted to spherical systems with isotropic v elocity dispersion ten- 
sors (e.g.. lTaffoni et alJl2003t lHavashi et alJl20fm The results of 
these simulations are also incorporated into semi-analytic mod- 
els which attempt to follow the tidal evolution of orbiting sub- 
haloes. However, the shapes and velocity dispersion tensors of 
haloes in cosmological simulations have been studied by many au- 
thors (e.g..lWarren et alJll992t ICole & Lacevlll996t | Thomas et alJ 
ll998tlColfn. Klypin & Kravtso\l200(tlBulloclj20021) and thev~ 
hibit a significant departure from spherical symmetry and isotropy. 
It is interesting to investigate if the response of a spherical subhalo 
to an imposed tidal field is the same as that of a triaxial subhalo. 

KMM compared the resilience to tidal forces of two identi- 
cal self-consistent NFW satellites on the same external tidal field 
and orbit, but constructed with different velocity dispersion ten- 
sors. Identifying the bound mass as function of time, these authors 
demonstrated th at the sate l lite with a r adially anisotropic velocity 
dispersion lOsipkov 1979; Merritt 1985) experiences much more 
efficient mass loss than its isotropic counterpart. This is due to the 
fact that particles on more radial orbits spend on average more time 
at larger radii and are therefore more easily stripped by the external 
field. Here we expand on the previous result and show that similar 
differences are obtained when one considers the tidal evolution of 
spherical and triaxial haloes as they orbit within a deeper poten- 
tial. We compare results of the 3:1 prolate and oblate haloes and 
the initial isotropic spherical halo used to create these models. The 
density profile of the spherical model is scaled by a factor of four 
in mass to preserve the same spherically averaged density profile 
as the prolate halo. Each halo is then placed on an identical or- 
bit within a static cored isothermal potential such that the ratio of 
host to satellite halo circular velocities is Vhost/fsat = snd 
the apocentric to pericentric radii is r apo /r pcr =4:1. The align- 
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Figure 8. Structural and kinematic properties of one subhalo in our highest-resolution cosmological run (R12) resolved with more than 2 X 10 5 particles. 
Left: The axial ratios b/a and c/a as a function of radius from the centre of the subhalo. Right: The anisotropy parameter /3 as a function of radius. In both 
cases, the thin lines correspond to measurements of the subhalo's properties just before entering the primary halo, whereas the thick lines correspond to its 
third pericentric passage. The subhalo becomes more spherical and isotropic as a result of the tides. Note that the radius of the subhalo is given in units of rt, 
where rt is its tidal radius after the third pericentre passage. 



ment of the axes of the triaxial satellite were varied with respect 
to the external deeper potential and the satellites always begin at 
the apocentre which is set equal to the virial radius. This allows the 
satellite-host system to be rescaled to any galactic or cluster sys- 
tem. The imposed potential is reasonably soft since we use a core 
radius equal to half the pericentric radius. Significantly higher rates 
of mass loss may be expected in a more cuspy potential, however 
here we aim to explore the slow tidal loss of mass rather than mass 
loss due to strong tidal shocks. In a forthcoming paper we will give 
an analytical model for the rate of mass loss as a function of host 
potential, orbit and subhalo structure and shape (Kazantzidis et al. , 
in preparation). 

In Figure [5] we perform a comparison of the tidal evolution 
of each of these three models by identifying the bound mass as a 
functio n of ti me using the publicly available group finder SKID 1 
fctadefeOOll) . The prolate halo experiences mass loss at a much 
higher rate than both the spherical and oblate ones. We also found 
that this result is not significantly sensitive to the alignment be- 
tween the major or minor axis of the satellite and the orbital plane 
with the maximum differences in the mass loss being of the order 
of 25%. In Figures |6| and we show the change in the shape and 
in the anisotropy parameter /3 of our models. We suspect that the 
prolate halo becomes more spherical due to the tidal removal of 
particles on radial orbits that are supporting the shape of the sys- 
tem. Future studies will analyse the orbital evolution of particles 
within the satellite to study these effects in detail. In addition, af- 
ter 5 Gyr of orbital evolution the velocity dispersion tensor of the 
prolate halo has significantly changed and now is nearly isotropic. 

It is interesting to study the evolution of shape and anisotropy 
for haloes in a CDM simulation that enter into a deeper potential. 
This is important when we want to model the orbits of stars in 
satellite galaxies or the effect of subhaloes on lense models. Here 
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we analyse subhaloes from a high resolution ACDM cluster sim- 
ulation ( Di emand et all2004l). The initial conditi ons are generated 
with the GRAFIC2 package (Bertschinger 2001). We start with a 
300 3 particle cubic grid with a comoving cube size of 300 Mpc 
(particle mass m p = 2.6 x lO lo /i _1 M0). We trace back and re- 
fine a cluster region by a factor of 12 in length and 1728 in mass, 
so that the mass resolution is m p = 1.5 x lO 7 ft _1 M0. The soft- 
ening length is comoving from the start of the simulation (z ~ 40) 
to z = 9. From 2 = 9 until present, we use a physical softening 
length of 1.05 x 10 _3 r v i r , paron t. At z = the refined cluster con- 
tains 1.4 x 10 7 particles (R12) and we select a relatively massive 
subhalo which is resolved with more than 2 x 10 5 particles. 

Figure|8|shows the axial ratios b/a and c/a (left panel) and the 
anisotropy parameter ft (right panel) as a function of radius for this 
object. The thin lines correspond to the subhalo just before entering 
its host, whereas the thick lines correspond to the third pericentre 
passage of the subhalo. Note that the radius of the subhalo is given 
in units of r t , where r t is its tidal radius after the third pericentre 
passage. Figure [8] confirms that subhaloes become more spherical 
and their velocity distribution more isotropic owing to tidal effects. 
We then compare the mean axial ratios of field haloes versus sub- 
haloes in the same simulation. On average we find that subhaloes 
have axial ratios that are 30% larger (more spherical) than field 
haloes. 



4 CONCLUSIONS 

We have presented a method for constructing cuspy axisymmetric 
and triaxial A-body haloes based on merging isotropic equilibrium 
spherical haloes with varying amounts of angular momentum. In 
particular, we found that radial mergers produce prolate systems, 
while mergers on circular orbits produce oblate systems. This tech- 
nique has the benefit that it is based on the way in which haloes 
obtain their triaxiality in cosmological simulations, therefore the 
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anisotropy distributions and angular momentum of haloes are well 
motivated. 

Mergers between similar equilibrium spherical haloes at high 
resolution show that the resulting halo has the same density pro- 
file, independent of the angular momentum of the merger. This has 
two implications: (i) the density profile of the triaxial halo can be 
set by the choice of the density profiles of the progenitor haloes 
and (ii) similar mass mergers can not dramatically re-arrange the 
central density structure of DM haloes, albeit oblate haloes (high 
angular momentum merger remnants) have concentrations up to 
a factor of two lower than prolate haloes (low angular momen- 
tum merger remnants). This may explain the entire scatter in the 
distribution of halo concentrations from cosmological simulations 
( Bullo ck et all200ll : lEke et aj200lklwechsler et all2002h . It may 
also be the case that the least concentrated haloes host the low sur- 
face brightness discs since these galaxies may form in high angular 
momentum oblate haloes (e.g., O'Nell et al. 1997). These galaxies 
generally require low value s of the concentration when fit to cuspy 
halo models (e.g. ,|McGaugh & de Blokll998| :lvan den Bosch et alJ 
2000: lde Blok et alJ200 it Is waters et all2003h . 

As an application we considered the tidal evolution of triaxial 
subhaloes orbiting within deeper potentials. Haloes with identical 
spherically averaged density profiles and on identical orbits evolve 
self-similarly with time. Spherical haloes with isotropic velocity 
dispersion tensors suffer significantly less mass loss than radially 
anisotropic prolate ones. This is likely due to the fact that particles 
on radial orbits are easily stripped which also results in subhaloes 
becoming more spherical. In general, the subhaloes become more 
spherical and their velocity distribution more isotropic owing to 
tidal effects. This result has been confirmed by investigating the 
response to tides of both isolated satellite haloes orbiting within a 
static host potential and substructure haloes in the time dependent 
cosmological tidal field. We find that subhaloes in a cosmological 
simulation of a cluster are on average 30% rounder than their field 
counterparts. Galaxy subhaloes should be even closer to spherical 
since they spend longer being reshaped by the host potential. 

Our more realistic modelling of DM haloes is important for 
studies of the weak and strong lensing statistics attempting to dis- 
tinguish between competing cosmological models, the structural 
evolution and mass loss from substructure, the formation of tidal 
streams or the sizes of satellite haloes in galaxies and clusters. 
These results may also be important to incorporate within semi- 
analytic models that attempt to model the distribution of satel- 
lites in DM haloes. CDM haloes are generally radially anisotropic 
therefore they lose mass and are disrupted more quickly than the 
isotropic systems that have been generically considered in previ- 
ous studies. This may explain why semi-analytic models predict far 
more substructures than found in high resolution numerical studies 
(e.g., IZentner & Bullockll2003l James Taylor private communica- 
tion). 
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